Equation of state of additive hard-disk fluid mixtures: 
A critical analysis of two recent proposals 
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A detailed analysis of two different theoretical equations of state for a binary mixture of additive 
hard disks [C. Barrio and J. R. Solana, Phys. Rev. E 63, 011201 (2001); A. Santos, S. B. Yuste and 
M. Lopez de Haro, Mol. Phys. 96, 1 (1999)], including their comparison with Monte Carlo results, 
is carried out. It is found that both proposals, which require the equation of state of the single 
component system as input, lead to comparable accuracy when the same input is used in both, but 
the one advocated by Santos et al. is simpler and complies with the exact limit in which the small 
disks are point particles. 

PACS numbers: 61.20.-p, 05.70.Ce, 64.10.+h, 51.30.+i 
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I. INTRODUCTION 

Despite being in principle a simpler system, hard-disk 
fluid mixtures have received much less attention in the 
literature than fluid mixtures of hard spheres. This may 
well be tied to the fact that up to now no analytical so- 
lution to the Percus-Yevick equation has been found for 
even dimensionality. In any case, what this has meant 
is that fewer results are available for fluid mixtures of 
hard disks than for hard-sphere mixtures. In particular, a 
very scarce number of proposals for the equation of state 
(EOS) of these mixtures has been made @, |, |, @, ||, al- 
though the trend seems to be reversing recently, and even 
fewer simulations have been performed to assess the value 
of such proposals. In a recent paper, Barrio and Solana 
H proposed an EOS for a binary mixture of additive 
hard disks. Such an equation reproduces the (known) 
exact second and third virial coefficients of the mixture 
and may be expressed in terms of the EOS of a single 
component system. They also performed Monte Carlo 
simulations and found that their recipe was very accurate 
provided an also very accurate EOS for the single com- 
ponent system (in their case it was the EOS proposed 
by Woodcock H) was taken as input. The comparison 
with other EOS for the mixture available in the literature 
indicated that their proposal does the best job with re- 
spect to the Monte Carlo data. Among these other EOS 
for the binary mixture considered in Ref. j5|, only the 
one introduced by Santos et al. a few years ago B also 
shares with Barrio and Solana's EOS the fact that it may 
be expressed in terms of the EOS for a single component 
system. The aim of the present paper is to present a 
detailed analysis of these two different equations of state 
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since the comparison made in Ref. || may be misleading 
in that it was not performed by taking the same EOS 
for the single component system in both proposals. A 
preliminary report of this work can be found in Ref. §. 

In order to carry out the analysis, the paper is orga- 
nized as follows. In Sec. II we recall the two different 
formulations for the EOS of a binary mixture of additive 
hard disks in a unified notation as well as provide the ex- 
plicit (approximate) expressions for the EOS of the single 
component system that will be used in the actual calcu- 
lations. This is followed in Sec. Ill by a discussion of the 
results and some concluding remarks. 



II. THE EQUATION OF STATE OF A BINARY 
MIXTURE OF ADDITIVE HARD DISKS 



Let us consider a binary mixture of additive hard disks 
of diameters <j\ and <72- The total number density is p, 
the mole fractions are x\ and xi = 1 — x\, and the pack- 
ing fraction is 77 = jp(a 2 ), where (a n ) = J2i x i a i- Let 
Z = p/pksT denote the compressibility factor, p be- 
ing the pressure, T the absolute temperature and ks the 
Boltzmann constant. Then, Barrio and Solana's EOS for 
a binary mixture of hard disks, Z^(rf), may be written 
in terms of a given EOS for a single component system, 
Z B (r)), as 

z* s (v) = 1 + + (i + tf [ZM 1] , (i) 

where £ = (a) 2 /(a 2 ) and /3 is adjusted as to repro- 
duce the exact third virial coefficient for the mixture -B3, 
namely 
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Here, 63 = 4(4/3 - VI/tt) is the reduced third virial 
coefficient for the single component system and B3 is 
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given by 



III. DISCUSSION 



B 3 = g {aiix 1 cr 1 + 3ai 2 a; 1 a;2cr 12 



(3) 



where 



ay = 7T + 2 (crf/cr?- - l) cos 1 (a l /2a lj ) 

-(44/a^l) 1/2 (l + af/24)^/2a^ (4) 

and <Jij = (oi + aj)/2. 

The EOS for the mixture, consistent with a given EOS 
for a single component system, introduced recently by 
Santos et al. reads 0| 
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We stress the fact that Eq. (0) is simpler than Eq. (|l|) 
[which must be complemented with Eqs. (0)— (Q)]. In ad- 
dition, the structure of Eq. (^J) is valid for any number of 
components, while Eq. (Q) requires the third virial coef- 
ficient, which is known exactly only for binary mixtures. 
In order to proceed with a quantitative analysis of Eqs. 
(|j) and (H), we have to specify Z s (rf). While many choices 
are available, we will restrict ourselves to the three fol- 
lowing EOS of the single component system: 

a) Woodcock's EOS [| 



1 - ^/??o 



+EM" 1 -4)(^r\ (6) 



n=2 



where rj = (\/3/6) 7r is the value of the crystalline 
close-packing and the b n (n = 2, . . . , 6) are the 
(known) reduced virial coefficients 

b) The Levin approximant of Erpcnbeck and Luban 



zf h (v) 
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71=0 
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where g„ = (-1)" \® j (1 - n/6) 5 b 6 /b e - n and 

Pn Em=0 ^7i+ 1 — mQm • 

c) The EOS proposed by Santos a/. @, H 



Vo 



(8) 



The EOS (|) and (@) are more complex than the EOS (|) 
in the sense that they require the exact knowledge of the 
first six virial coefficients, while Eq. (^) is constructed by 
using the first two virial coefficients only and enforcinga 
pole at rj = t]q. Despite its simplicity, however, Eq. (g) 
does a remarkably good job when compared with simula- 
tion data, although it is of course less accurate than the 
more sophisticated EOS (§) and (R) JlO| . 



In Table Q, we show the results of Eqs. (|l|) and (||) 
when Woodcock's EOS 0, Eq. (^), for the single com- 
ponent system is used as input in both equations, as well 
as the available MC data H. As seen in Table |, it is 
fair to say that both recipes are of comparable accuracy 
with respect to the Monte Carlo results, their difference 
being generally smaller than the error bars of the sim- 
ulation data, although Z^(rf) performs slightly better 
than Z^ YH (r/). This may be fortuitous since if one takes 
for ZJrf) in Eqs. (|l|) and (||) the Levin approximant ||, 
Eq. ([7]), (which is known to give the most accurate ap- 
proximation to the single component compressibility fac- 
tor |9|, 10 ) the apparent (slight) superiority of Z^(n) 
2;er there. For instance, the theoretical values 
corresponding to the packing fraction 77 = 0.6 



is no long 
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are increased by about 0.03 when the Levin approximant 
rather than Woodcock's EOS is used as input, so that 
in this case the accuracy of (77) is generally slightly 
better than that of Z^(rf). This is further illustrated 
in Figs. [l|-|J where in order to enhance the differences 
we have plotted Z m (r}) — Z^ (77) versus 77, Z^ (77) 
indicating the use of Eq. |l]) taking for Z s (rj) the EOS 
by Woodcock. The figures show that, in general, the dif- 
ferences between the EOS (|l|) and the EOS (||), taking 
of course the same Z s as input, are smaller than or of 
the order of the error bars of the simulation data. As 
expected, a better agreement with the simulation data is 
obtained when either of the more accurate EOS @ or 
(0) is used as input instead of the much simpler EOS 
It is interesting to remark that the best agreement at the 
two largest densities, 77 = 0.55 and 77 = 0.6, corresponds 
to the use of the Levin approximant for Z s . 

Let us try to understand why both EOS for the mixture 
give practically equivalent results when the same input 
is used in both. First, it may be shown that ^ yh (t?), 
while not reproducing the exact third virial coefficient 
S3 , yields a very good estimate of it [O , namely -B3 ~ 



[1 + (6a (tt/4) 2 
into Eq. (g), we get 

/3~ 



a ) . If we replace that estimate 
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By using this estimate in Barrio and Solana's EOS, 
Eq. (nl), we have 



^ s (r7)-^(r7)_.(l-OA(r7), 
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where 



AM s - 2 



>-ll-I 



[z s (v) - 1] - 



(10) 
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According to the approximation involved in Eq, 
the difference Z^(rf) — Z^ i{ (ff) is small if the asym- 
metry of the mixture is small (£ < 1) and/or A(?7) is 
small. The function A(rj) is plotted in Fig. |] for the 
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TABLE I: Compressibility factor for different binary mixtures of hard disks as obtained from Monte Carlo simulations, from 
Eq. ([!]), and from Eq. In the two latter, Woodcock's equation of state for the single component system is used. 



Xi = 0.25 



xi = 0.50 



xi = 0.75 



C2/01 


V 


MC a 


2/3 


0.20 


1.559(6) 




0.30 


2.036(8) 




0.40 


2.79(1) 




0.45 


3.31(1) 




0.50 


4.02(2) 




0.55 


4.98(2) 




0.60 


6.31(2) 



Eq. (1) Eq. (5) 



MC a 



Eq. (1) Eq. (5) 



MC a 



Eq. (1) Eq. (5) 



1.559 
2.040 
2.79 
3.32 
4.02 
4.97 
6.29 



1.559 
2.040 
2.79 
3.32 
4.02 
4.97 
6.28 



1.561(5) 
2.043(8) 
2.79(1) 
3.31(2) 
4.02(1) 
4.98(2) 
6.30(1) 



1.558 
2.039 
2.78 
3.32 
4.02 
4.96 
6.28 



1.558 
2.039 
2.78 
3.32 
4.02 
4.96 
6.27 



1.565(4) 
2.051(8) 
2.80(1) 
3.33(1) 
4.04(2) 
5.03(2) 
6.36(3) 



1.563 
2.048 
2.80 
3.34 
4.05 
5.00 
6.33 



1.563 
2.048 
2.80 
3.34 
4.05 
5.00 
6.33 



1/2 


0.20 


1.534(6) 


1.536 


1.536 


1.540(6) 


1.538 


1.538 


1.556(7) 


1.552 


1.552 




0.30 


1.998(7) 


1.995 


1.995 


2.008(7) 


2.000 


2.000 


2.039(8) 


2.027 


2.026 




0.40 


2.72(1) 


2.70 


2.70 


2.71(1) 


2.71 


2.71 


2.77(1) 


2.76 


2.76 




0.45 


3.20(2) 


3.21 


3.21 


3.22(2) 


3.22 


3.22 


3.29(2) 


3.29 


3.29 




0.50 


3.88(1) 


3.88 


3.87 


3.90(2) 


3.89 


3.89 


3.98(2) 


3.98 


3.98 




0.55 


4.79(2) 


4.77 


4.76 


4.81(2) 


4.80 


4.78 


4.93(2) 


4.91 


4.90 




0.60 


6.03(3) 


6.02 


6.00 


6.04(2) 


6.05 


6.03 


6.22(3) 


6.21 


6.20 


1/3 


0.20 


1.491(6) 


1.490 


1.490 


1.510(8) 


1.506 


1.506 


1.538(9) 


1.536 


1.536 




0.30 


1.907(8) 


1.905 


1.904 


1.940(8) 


1.937 


1.936 


2.004(9) 


1.996 


1.995 




0.40 


2.55(1) 


2.54 


2.54 


2.59(1) 


2.60 


2.60 


2.71(1) 


2.71 


2.70 




0.45 


2.99(2) 


3.00 


2.99 


3.07(1) 


3.07 


3.06 


3.20(2) 


3.21 


3.21 




0.50 


3.60(2) 


3.59 


3.57 


3.69(2) 


3.69 


3.68 


3.89(2) 


3.88 


3.87 




0.55 


4.39(3) 


4.39 


4.36 


4.52(2) 


4.52 


4.50 


4.79(2) 


4.78 


4.76 




0.60 




5.49 


5.44 


5.66(6) 


5.68 


5.64 


6.06(1) 


6.03 


6.00 



a Ref. Q 



cases where Z s (rj) is given by Woodcock's EOS, by the 
Levin approximant, and by the SHY EOS. In all in- 
stances it is practically zero up to r\ sa 0.2 but then 
it grows rapidly. The most disparate mixture consid- 
ered in Barrio and Solana's simulations corresponds to 
a; 1 = 0.25, a = <Jijo\ = 1/3, which yields £ = 0.75. 
This explains the fact that Z* s (r)) - Z^(r)) < 0.05 
in the simulated cases. It should be noted however 
that Eq. ( JlCj) tends to overestimate the actual differ- 
ence Z^(rj) — Z^ yh (t7) (for instance, in the case where 
x\ = 0.25, a = 1/3, 77 = 0.6 and Z s (j]) is given by Wood- 
cock's EOS, this difference is 0.05, while the prediction of 
Eq. ( [l0| ) yields 0.07), so that its main purpose is to illus- 
trate the fact that both EOS yield practically equivalent 
results for not very asymmetric mixtures. On the other 
hand, more important differences can be expected for 
disparate mixtures, especially in the case of large densi- 
ties. At a given density and a given diameter ratio a < 1 
the smallest value of the parameter £ corresponds to a 
mole fraction x\ = a/(l + a) for the large disks, namely 
£ = 4a/ '(1 + a) 2 . Thus, £ <C 1 if a <C 1 and, according 
to Eq. ©, Z* S ( V ) - Z£™(r,) ~ A(tj). 

Let us consider now the limit in which the small disks 
become point particles (a — > 0) and occupy a negligible 



fraction of the total area. In that case, the compressibil- 
ity factor of the mixture must reduce to [|| p"3[ 



X-2 



n 



xxZ s {rj). 



(12) 



The first term represents the (ideal gas) partial pressure 
due to the point particles in the available area (i.e., the 
total area minus the area occupied by the large disks), 
while the second term represents the partial pressure as- 
sociated with the large disks. In the limit a — > with X\ 
finite (or, more generally, for a <C x\), we have £ — > x\ 
and j3 — ► (1 — ba/2)x2/(^ + %i) [note that in this limit 
the approximation (||) becomes correct], so that 



7SYH 



(V) 



X-2 



I + X1+I2 1 



■xiZ 8 (r)), 



h 
2 



(13) 



[ZM - 1] 



(14) 

Therefore, while Eq. (|5|) is consistent with the exact prop- 
ertyjjl^), Eq. ([!]) violates it. In fact, the right-hand side 
of fllQ|), with £ = xi, gives the deviation of Barrio and 
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FIG. 1: Plot of the difference 8Z(j)) = Z m (n) - Z^ S(W) (77) 
for a size ratio aifax = 2/3 and for x\ — 0.25 (top panel), 
x\ = 0.50 (middle panel), and xi = 0.75 (bottom panel). The 
filled circles are Monte Carlo data [^|, the dashed lines refer 
to the proposal (|l|), and the solid lines refer to the proposal 
(^). The EOS of the single component fluid, ZJri), used as 
input are Eq. Q (lines without symbols), Eq. (Q) (lines with 
squares), and Eq. (0) (lines with triangles). 



Solaria's EOS from the exact compressibility factor in the 
special case a — > 0. 

In summary, in this report we have performed a de- 
tailed comparison of the EOS proposed by Barrio and 
Solana |3| and the one introduced by Santos et al. ||. 
We find that both proposals lead to comparable accuracy 
when the same EOS for the single component system is 
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FIG. 2: Same as in Fig. [I], but for a size ratio <J%joi = 1/2. 



used and confirm that the more accurate the Z s {rf) the 
more accurate the resulting compressibility factor for the 
binary mixture. In favor of Santos et aL's proposal, apart 
from its simpler form which also yields a very reasonable 
estimate of the known third virial coefficient, is the fact 
that it is readily extendible to the multicomponcnt case 
(including polydisperse mixtures) and complies with the 
exact limit in which the small disks are point particles, 
while Z^(rj) does not share these assets. 
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FIG. 4: Plot of A(»7), Eq. (0), by assuming Woodcock's EOS 
(dashed line), the Levin approximant (solid line), and the 
SHY EOS (dotted line) for the pure fluid. 



